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', Abstract 

^ . We discuss the potential of future sub-10 GeV threshold imaging atmospheric 

^T) I Cherenkov telescope arrays for exploring the physics of rotation powered pulsars 

' and their interactions with the ambient medium through relativistic winds and 

. termination shocks. One such telescope is the high-altitude concept called "5@5" 

^ I recently suggested by Aharonian et al. (2001). 5@5, with its enormous detection area 

' exceeding 10^ at the threshold energy of about 5 GeV, combines two distinct 

• features of the current satellite-borne (large photon fluxes at GeV energies) and 

^p^i ground-based (large detection areas at TeV energies) gamma-ray astronomies. Such 

O I an instrument would allow comprehensive studies of temporal and spectral char- 

> ' acteristics of 7-ray pulsars in the crucial 5 to 30 GeV energy interval. An equally 

. important topic in the program of pulsar studies by 5@5 would be the search for GeV 

^ I 7-radiation from other radio pulsars at a few mVela level. And finally, the searches 

for pulsed radiation components in the spectra of a large fraction of unidentified 
EGRET sources (suspected to be pulsars) without invoking information from lower 
I (radio, optical. X-ray) frequency domains, seems to be another important issue, 

because the periodic signals at lower energies could be significantly suppressed in 
many cases. The detection rate of 7-rays from "standard" EGRET sources by 5@5 
is expected to exceed several events per one second. This should provide an ad- 
equate photon statistics for the search for periodic signals at the flux level of 3 
mVela within the observation time of 3 h or so (a time resolution below which any 
change of a signal's phase can be ignored). The spectral coverage by 5@5 and its 
flux sensitivity are nicely suited for studying other aspects of pulsar physics and 
astrophysics, in particular for detecting unshocked relativistic pulsar winds, as well 
as for quantifying characteristics of pulsar driven synchrotron nebulae through the 
inverse Compton radiation at energies between several GeV and several 100 GeV. 
The Vela pulsar, the brightest 7-ray source on the sky, is an ideal laboratory for 
practical realization of these unique observational possibilities. 
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1 Introduction 

The number of cataloged radio pulsars - single neutron stars powered by fast 
rotation - exceeds 1200 (Lorimer 2001). Only seven of them (with two other 
possible candidates) are reported by the EGRET team as GeV 7-ray sources 
(Thompson 1999, Kaspi ct al. 2000) . Actually, this is not a big surprise, and 
can be explained by modest flux sensitivities of 7-ray instruments. For ex- 
ample, the minimum detectable flux by EGRET of about 10~^^ erg/cm^s is 5 
orders of magnitude larger than the sensitivity achieved by radio telescopes 
at GHz frequencies. This striking difference in sensitivities is compensated, 
to a certain extent, by much larger energy fluxes in 7-rays. While the pul- 
sars emit in the radio band a tiny < 10~^ fraction of their rotational energy 
(Manchester & Taylor 1977), the 7-ray luminosities of some of the EGRET 
pulsars above 100 MeV exceed one per cent of their spin-down luminosities. It 
is expected that the Gamma-ray Large Area Space telescope (GLAST), the 
future major satellite-based high energy 7-ray detector with flux sensitivity as 
good as 10"^^ erg/cm^s (Gehrels & Michelson 1999), will increase the number 
of radio pulsars seen also in 7-rays by at least an order of magnitude (Thomp- 
son 2001). Moreover, the radio pulsars, which are observable only when their 
radio beams cross the Earth orbit, are only a part of a larger source popu- 
lation called Rotation Powered Pulsars (RPPs, for short). The 7-ray beams 
are believed to be signiflcantly wider than the radio beams and they should 
not necessarily overlap with each other (Yadigaroglu & Romani 1995). If so, 
the RPPs would have improved chances, in principle, to be detected, in high 
energy 7-rays rather than at radio wavelengths. This could be the case of some 
of the unidentifled EGRET sources, a large fraction of which is believed to 
be radio-quiet pulsars (see e.g. Grenier 2000). The unambiguous and straight- 
forward proof of this hypothesis would be the discovery of 7-ray pulsations 
from these objects. In this regard, it is crucial to search for the pulsed 7-ray 
components without relying on observations at other energy bands. Although 
GLAST will be able to perform effective searches for strong unidentifled 7-ray 
sources (Thompson 2001), its potential in this regard will be still limited, at 
least for marginally detected weak sources which most likely will appear in the 
GLAST source catalog (future "unidentifled GLAST sources"). An effective 
realization of this important approach which relies only on 7-ray astronom- 
ical observations requires large photon statistics, and therefore much larger 

^ E-mail: Felix.Aharonian@mpi-hd.mpg.de 
^ E-mail: bogoval@axpk40.mephi.ru 



2 



detector areas. The current and forthcoming atmospheric imaging Cherenkov 
telescopes do provide huge (as large as 10^ m^) detection areas. But most 
of these telescopes operate in the energy region > 100 GeV which lies, most 
probably, well beyond the cutoffs expected in the pulsar spectra. In contrast, 
the recently suggested concept "5@5" by Aharonian et al. (2001) -5 GeV en- 
ergy threshold array of imaging atmospheric Cherenkov telescopes at 5 km 
altitude- which combines the most important features of the current space- 
based (large fluxes at GeV energies) and ground-based (large detection areas 
at TeV energies) astronomies, could serve as an ideal tool for pulsar studies. 
Remarkably, the analysis of the EGRET pulsars shows that with age of pul- 
sars the fraction of the spin-down luminosity converted into 7-rays increases, 
and the peak in the spectral energy distribution vSp shifts towards 10 GeV 
(Nel et al. 1996). If this tendency extends to fainter pulsars, the chances for 
detection of many pulsars by 5@5 would be increased dramatically. 

The potential of the 5@5 concept is not limited to the discovery of faint 7- 
ray pulsars. In addition, 5@5 can provide detailed spectroscopy in one of the 
key or, perhaps even the most informative, energy regions of above several 
GeV for pulsar physics. It is worth noting that the EGRET data indicate that 
some important changes might take place at multi-GeV energies. In partic- 
ular, the light curves of all EGRET pulsars at energies above 5 GeV seem 
to be different than at lower energies (Thompson 2001). Detailed temporal 
and spectroscopic measurements in this transition region by 5@5 will hope- 
fully remove many uncertainties we presently face in the physics of pulsar 
magnetospheres. Accurate measurements of the energy spectrum of strong 
EGRET pulsars in the cutoff region expected between 5 and 20 GeV would 
allow, in particular, to distinguish between two currently popular models of 
high energy 7-radiation — polar cap (Daugherty & Harding 1982, 1996; Usov 
& Melrose 1995) and outer gap (Cheng et al. 1986; Romani 1996; Hirotani 
& Shibata 1999) models. These models generally predict different fractions 
in the population of radio-quiet versus radio-loud 7-ray pulsars, as well as 
different spectral characteristics, especially at the multi-GeV energy range. 
These differences are quite significant, even with large uncertainties in model 
parameters (Harding 2001). Thus 5@5 should allow distinction among these 
models, or perhaps even to challenge most scenarios. Remarkably, this inter- 
mediate energy region can be effectively studied by both GLAST and 5@5 — 
a combination of which would provide high reliability of results. On one hand 
GLAST will extend the studies to low energies, down to ~ 20 MeV, while 5@5 
can perfectly cover the energy region well above 10 GeV where GLAST most 
probably would run out of photons. Presently several ground-based telescopes 
such as MAGIC, CELESTE and STACEE are pushing their energy thresholds 
below 100 GeV (e.g. Buckley et al. 2001; Krennrich 2001), however none of 
these is optimized for energies in the crucial 10 GeV range. 

The flux sensitivity, angular resolution, and the spectral coverage of 5@5 per- 
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fectly match other important aspects of the pulsar physics and astrophysics. 
In particular, the search for inverse Compton 7-radiation from unshocked ul- 
trarelativistic winds caused by illumination of the wind by the IR to X-ray 
radiation emitted by the pulsar/neutron star itself (Bogovalov & Aharonian 
2000) or, in the case of binary pulsar systems, with illumination provided 
by the companion star (Ball & Kirk 2000), would be an important topic in 
the program of pulsar studies by 5@5. Because this radiation is produced in 
deep Klein-Nishina regime, it has a rather specific shape - unusually narrow 
spectra peaked at ~ F^meC^, where is the wind Lorentz factor. The 
detection of this radiation would be the first observational test for existence 
of ultrarclativistic pulsar winds, and would give unique information about the 
Lorentz factor and the site of creation of the kinetic energy dominated (KED) 
wind. The large dynamical energy range of 5@5, extending from several GeV 
to several hundred GeV, would allow such a search for pulsar winds in the 
most relevant range of Lorentz factors between 10^ to 10^. 



Another, more traditional, objective for 5@5 would be the detection of inverse 
Compton 7-rays from synchrotron nebulae surrounding pulsars - the regions 
powered by the wind termination shocks (plerions). The TeV radiation de- 
tected from the Crab Nebula by many groups (for review see e.g. Catancse 
& Weekes 1999) agrees reasonably well with model calculations (De Jager & 
Harding 1992; Dc Jager et al. 1996; Atoyan & Aharonian 1996; Aharonian & 
Atoyan 1998; Hillas et al. 1998) performed within the MHD model of Ken- 
nel & Coroniti (1984). According to this model, the wind is terminated by a 
standing reverse shock at a distance of about r = 0.1 pc, the shock in turn ac- 
celerates electrons up to energies exceeding 10^^ eV and randomizes their pitch 
angles. The inverse Compton 7-radiation at GeV/TeV energies, combined with 
synchrotron optical and X-ray emission, contains information about the rel- 
ativistic electrons and the nebular magnetic fields in the downstream region 
of the shock. But the expected 7-ray fluxes from other pulsar-driven nebu- 
lae, based on the observed synchrotron X-ray fluxes as well as on the model 
assumptions concerning the nebular magnetic flelds, contain signiflcant un- 
certainties. Nevertheless, we believe that the future low-energy lACT arrays 
should be able to reveal and allow study of the spectral and spatial character- 
istics of this component of radiation from the synchrotron nebula surrounding 
the Vela pulsar, and hopefully also from some other pulsars with "spin-down 
fluxes" E/A-n-d"^ > 10^^ erg/s kpc^ (provided that a noticeable part of the spin- 
down luminosity can be transformed eventually into the shock-accelerated TeV 
electrons) . 
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2 Performance of 5@5 



The concept of 5@5 is basically motivated by the necessity of sensitive "gamma- 
ray timing explorer" for detailed studies of the characteristics of transient 
GeV 7-ray phenomena (Aharonian et al. 2001). In this regard, 5@5 and 
GLAST present two complementary approaches. While GLAST with its al- 
most 27r steradian field of view will allow an effective monitoring of large 
number of persistent (quasi-stable) 7-ray sources, 5@5 can provide a deeper 
probe of individual highly variable or transient sources. For example, 5@5 
can detect any multi-GeV 7-ray flare with an apparent luminosity^ of about 
10^^((i/l Gpc)^erg/s at < 10 GeV {d is the distance to the source normal- 
ized to 1 Gpc), lasting only a few hours. Therefore the prime objective of 
instruments like 5@5 will be studies of high energy transient phenomena like 
GeV/TeV flares of blazars (and perhaps also microquasars) as well as solitary 
gamma-ray events like GRBs (Aharonian et al. 2001). The second important 
objective of 5@5 concerns the detailed 7-ray spectroscopy of persistent 7-ray 
sources in the energy region between 10 GeV and 100 GeV. 

In the original paper on the concept "5@5" (Aharonian et al. 2001), the fol- 
lowing conflguration for the telescope array has been proposed. The system 
consisting of 5 lACTs is to be installed at the altitude H = 5 km a.s.l., close to 
the ALMA (the Atacama Large Millimeter Array) site in Northern Chile. Four 
telescopes are located at the corners, and one in the center of a square with 
a linear size d — 100 m; each telescope has an approximately 20-30 m diame- 
ter optical reflector, and is equipped with a multi-channel conventional PMT 
based camera with a pixel size 0.12° and an effective held of view FoV= 3.2°. 



2.1 Detection area, angular and energy resolutions 

To understand the potential of 5@5 we need two basic parameters - the ef- 
fective detection area Ae{{{E) and the angular resolution (f){E) as functions of 
photon energy, approximated in the range from 1.5 GeV to 100 GeV in the 
following convenient analytical forms: 

Aff = 8.5E^-^[1 + {E/5 GeV)^-']-^ m' , (1) 



^ In the case of sources with strong beaming (Hke pulsars) within a solid angle 
Q < 1 ster, or sources moving relativistically towards the observer (like blazars) 
with Doppler factor (5j 3> 1, the intrinsic 7-ray luminosities are smaller than the 
apparent luminosities by factors U/4ir and respectively. 
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and 



= 0.42(^/5 GeV)"°-Megree, (2) 

where E is the energy of the primary 7-ray photon or electron (Aharonian et 
ah 2001). The detection area has a strong energy-dependence below 10 GeV, 
Aes{E) (X but at higher energies the detection area increases slowly with 
energy, Aef[{E) oc E^-^. Close to 1 TeV the detection area actually becomes 
constant, in essence because of the limited field of the view of the camera. 
The parameter 0, which corresponds to the half-angle of the cone centered 
on the source and containing 67% of events at energy E, defines the energy 
dependence of the Point Spread Function (PSF) of the telescope. 

The energy resolution of the instrument depends on the energy; at energies 
above 10 GeV it could be close to 15-20 per cent. In this paper for simplicity 
we assume an average estimate for the energy resolution, namely a constant 
value at the level of 20 per cent throughout the entire energy interval. 

The above numbers characterizing the basic parameters of 5@5 are based on 
the preliminary Monte Carlo calculations (Aharonian et al. 2001) and need 
to be confirmed by independent simulations, especially in the sensitive energy 
region of primary 7-rays below 10-20 GeV where the atmospheric cascades 
consist of only a few generations of secondary particles. Further note that the 
assumed configuration of telescopes is based on our "best guess". The opti- 
mal configuration of both the individual telescopes (the area of the dish and 
the field of view and pixel size of the imager) and the system (the number 
of units and average distance between the telescopes) is a subject for future 
comprehensive Monte Carlo studies. One cannot exclude that the new perfor- 
mance studies may result in some deviations from the parameters used in this 
paper. However, we believe that the possible uncertainties in the flux sensi- 
tivities should not exceed a factor of 2 or 3. Therefore, the main conclusions 
of this paper concerning the potential of the "5@5" concept for pulsar studies 
are rather independent of specific configurations of future sub- 10 GeV lACT 
arrays. 

2.2 Background conditions 

At energies below 30 GeV the background detected by the system of lACTs 
operating in the stereoscopic mode is dominated by electromagnetic showers 

produced by cosmic ray electrons. The dominance of the electronic compo- 
nent is explained by the combination of several factors, in particular (i) by 
the large, up to a factor of 5 to 10 difference between the energies of electrons 
and protons producing the same amount of Cherenkov light, (ii) by the high 
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Fig. 1. The energy distribution of background events accumulated during obser- 
vation time Tabs = W^s within the sohd angle AQ = 10~^sr, contributed from 
the following channels: (1) by cosmic-ray electrons (dotted curves curves marked 
as "0", "10", and "20" correspond to geomagnetic cutoffs at 0, 10, and 20 GeV, 
respectively); (2) from the diffuse galactic 7-ray background in the direction of the 
galactic center (GD(C)) and anticenter (GD(A)); (3) from the diffuse extragalactic 
7-ray background (EGB); from (4) a "typical" molecular cloud (GMC), and (5) a 
"typical" supernova remnant (SNR) (see the text for details). 

altitude of observations, (iii) by compact (compared with hadronic showers) 
Cherenkov images of electromagnetic showers, (iv) by significant increase (ap- 
proximately oc E~^'^) of the electron-to-proton ratio of cosmic rays down to 
E GeV. This offers a unique method for the continuous {on-line) calibra- 
tion and control of the detector characteristics using the cosmic ray electrons 
(Aharonian et al. 2001). The differential spectrum of cosmic ray electrons in 
the region from several GeV to 1 TeV can be approximated as: 

= 1.36 X 10-3 E-\l + (E/5GeV)2-2]-i cm-^s'^sr-^GeV-^ (3) 

ail 



The experimental data used for this approximation contain large, up to factor 
of 2 uncertainties, not only at very high energies, but also in the most criti- 
cal energy region around 10 GeV. Moreover, in this energy region the flux of 
electrons entering the Earth's atmosphere depends on the geomagnetic field, 
and therefore on the detector location. For the ALMA site with coordinates 
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23° S and 67° W, the geomagnetic cutoff of the electrons in the vertical di- 
rection is about 13 GeV (Quenby & Wenk 1962). This value agrees with the 
analytical prediction for the dipole magnetic field with the equator shifted to 
the south to 10° at this longitude (Adams et al. 1981). The estimates in the 
dipole approximation give geomagnetic cutoffs in the electron spectra close 
to the horizon in the eastern and western directions at 9.4 GeV and 33 GeV, 
respectively. 

In Fig. 1 we show the detection rates of primary electrons by 5@5 (dotted 
curves) calculated for 3 different assumptions concerning the geomagnetic cut- 
off - 0, 10 GeV, and 20 GeV. For the ALMA site the geomagnetic cutoff always 
exceeds 10 GeV; for a broad range of optimal (for Cherenkov observations) 
zenith angles < 40°, the geomagnetic cutoff varies between 10 and 20 GeV. 
Since the electron spectrum below the cutoff has a very sharp (almost trun- 
cated) form, the cosmic ray background in this energy region is dominated 
by protons of energy above the geomagnetic cutoff. It should be noted, how- 
ever, that the geomagnetic cutoff is relevant only to the primary electrons, 
but not to the secondary ones produced in the atmosphere by effects of high 
energy protons. Therefore we should expect electron fluxes extending well 
below the geomagnetic cutoff. The recent measurements of this component 
of electrons at energies below 10 GeV by the Alpha Magnetic Spectrometer 
(AMS) indeed revealed non-negligible fluxes comparable with the flux of pri- 
mary electrons (Alcaraz ct al. 2000). Thus this component cannot be a priori 
neglected in estimations of potential contributors to the background for space 
and ground based gamma-ray instruments. However, for 5@5 type of instru- 
ments the contribution of this component fortunately appears to be quite 
modest, namely less than the contribution from hadronic showers induced by 
high energy (above the geomagnetic cutoff) primary protons, but accepted by 
the telescopes at both the hardware (trigger) and software (image identifica- 
tion) levels as low (sub geomagnetic cutoff) energy electromagnetic showers. 
The reason is that the trapped electrons spend most of their time in the vol- 
ume where the AMS did measurements (e.g. Lipari 2002), therefore the flux of 
secondary electrons close to the Earth's atmosphere is estimated significantly 
(by an order of magnitude) lower than the fluxes observed at the height of the 
AMS orbit (Lipari, private communication). 

At an observation altitude of 5 km a.s.l. the 30-100 GeV protons produce 
approximately the same amount of Cherenkov light than the 7-rays and elec- 
trons with energy 10 GeV. Therefore, the cosmic ray protons with energy 
above the geomagnetic cutoff initiate showers a non-negligible part of which 
could be accepted as low energy electromagnetic showers. The detailed study 
of such "small-proton" events requires detailed Monte Carlo simulations. Our 
preliminary study shows that the integral rate of such events does not exceed 
10 per cent of the detection rate of cosmic ray electrons. Moreover, the further 
analysis of the shapes of these events should allow an additional (by a factor 
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of 3 or so) suppression of the proton background. Therefore in Fig. 1 we as- 
sume that below the rigidity cutoff the cosmic ray background rate does not 
disappear but stays at the level of several per cent of the "nominal" (i.e. with- 
out rigidity cutoff) electron detection rate. Such an approximation should be 
confirmed by future detailed studies of the misidentified low-energy hadronic 
showers. Therefore, the detection rates shown in Fig. 1 by the curve with- 
out geomagnetic cutoff and by two curves with geomagnetic cutoff at 10 and 
20 GeV should be considered as conservative upper and lower limits for the 
cosmic ray background. 

The isotropic diffuse extragalactic background (EGB) measured by EGRET 
in the energy interval 5-300 GeV is approximated as (Sreekumar et al. 1998) 

'^'^''wo^^^ ^ ■ GeV)-2-^5 ph/cm^s sr GeV . (4) 



The fluxes of diffuse galactic 7-radiation, produced by interactions of cosmic 
rays with the interstellar medium, and partly contributed by unresolved galac- 
tic sources, strongly depend on the galactic coordinates. The largest flux of 
this component arrives from the direction of the galactic center. For example, 
the average flux from the region I — 0° and 2° < |6| < 6° can be approximated 
as (Hunter et al. 1997) 

~ 3 . 10-6(^/5 GeV)-'-^ ph/cm's sr GeV . (5) 

ail 



The diffuse background from the region / = 180°, 6° < b < 10°, i.e. in the 
direction of the Anticenter, is approximated as (Hunter et al. 1997) 

= 1.45 • 10-'(E/5 GeV) -2-^3 ph/cm^s sr GeV . (6) 



In the case of location of an weak 7-ray pulsar close to (or behind) a rel- 
atively young supernova remnant (potential sites of cosmic ray accelerators 
and, therefore, potential high energy 7-ray emitters), or giant molecular clouds 
(radiating 7-rays due to interactions with the "sea" of galactic cosmic rays) 
additional contamination from these extended sources (with angular radius 
up to 1° or so) should be taken into account (see e.g. Aharonian 1995). For a 
"typical" 7-ray emitting SNR, we use the following approximations 

Jsnr{E) = 4 • 10-^(^/5 GeV) -2 ph/cm^s GeV (7) 



where Ai = W^5o^i/'^kpc) = VFp/10^°erg is the total energy of accelerated 
protons contained in the remnant, Ui = n/lcm~^ is the density of ambient 
gas, and dkpc = d/1 kpc is the distance to the source. 
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For a "typical" GMC we adopt the flux 

Jmc{E) = 5.6 • 10-^°(^/5 GeV)-^-^^ A2 ph/cm^s GeV , (8) 

where A2 = M^/dl^^- M5 = M/IO^M© is the the mass of the molecular cloud. 

The detection rates of cosmic ray electrons, as well as of the diffuse galactic 
and extragalactic 7-ray backgrounds within Af2 = 10~^ sr solid angle (opening 
angle ~ 1°) are calculated with 

R{E) = j J{E')Aeff{E')AnM{E - E')dE' , (9) 

where M(E — E') = exp((E — E')'^ /2a'^) / (\^2TTa) is the energy spread function 
for the Gaussian distribution with a — E/5 (energy resolution 20%). 

For the SNR and GMC type of sources with "typical" values for the scaling 
parameters Ai = 0.1 and A2 = 1, respectively (Aharonian 1995), and assuming 
~ 1° angular radius (solid angle Afl = 10~^), 

R{E) = J J(E')Aeff(E')M{E - E')dE' . (10) 

In Fig. 1 we show the differential energy distributions of detected background 
events. The corresponding total number detected counts accumulated during 
the exposure time Tobs = 10^ s (typical observational time of a single source 
by Cherenkov telescopes during one night) in the energy range 3 to 300 GeV 
are: (i) from cosmic rays - 8.1 ■ 10^ and 2.8 • 10^ for geomagnetic cutoffs at 10 
GeV and 20 GeV, respectively; (ii) 5.6 • 10^ and 3.2 • 10^ counts for the diffuse 
7-ray galactic emission from directions to the Galactic Center and Anticenter, 
respectively; (iii) 1.2 ■ 10"^ counts for the diffuse extragalactic background; 
(iv) 1.3 ■ 10^ and 9.4 ■ lO''^ counts from "typical" SNRs and GMCs. These 
results show that except for the most "contaminated" directions towards the 
Galactic Center, the cosmic ray induced background always dominates over 
the contributions from the extragalactic and galactic 7-ray backgrounds. 



3 High energy gamma-rays from RPPs 

High energy 7-rays from RPPs are produced in 3 physically distinguished 
regions through several radiation mechanisms. These sites are schematically 
shown in Fig. 2. The pulsed 7-rays are generated in the pulsar magnetosphere 
limited by the light cylinder. The relativistic wind ejected by the pulsar can 
itself produce 7-rays as well as to excite the surrounding nebula which emit 



10 



Radiation from a Pulsar-wind-nebula complex 




Pulsar 
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6 e\ Shock front 



Synchrotron 
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y: MeV/GeV/TeV 



Interstellar medium 




Fig. 2. Sketch of the sites and radiation mechanisms of nonthermal emission asso- 
ciated with rotation powered pulsars: (i) the region within the light cylinder where 
the magnetospheric pulsed radiation from radio to 7-rays is produced; (ii) the part 
of the wind of cold relativistic plasma close to the light cylinder which effectively 
emits GeV and TeV 7-rays through the IC mechanism; (iii) the surrounding syn- 
chrotron nebula (plerion) which emits broad-band electromagnetic radiation from 
radio to multi-TeV 7-rays through the synchrotron and IC channels. 

electromagnetic emission in a wide electromagnetic band. Gamma-rays from 
these regions carry crucial information about the pulsar and its interaction 
with the surrounding medium. 



The brightness temperature of the pulsed radio emission is of the order of 
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1023 _ 1030 K (Manchester & Taylor 1977) implies that this emission is gen- 
erated coherently by beams of relativistic electrons. In the IR/optical and 
X-ray bands the emission is not coherent; it may have thermal or nonthermal 
origin. The specific mechanism(s) responsible for the radio (e.g. Lyubarsky 
1995) as well as IR/optical and X-ray (Pacini 1971, Becker & Triimper 1997, 
Yadigaroglu & Romani 1995, Harding & Zhang 2001) remain highly uncertain 
even for the most prominent objects - the Crab and Vela pulsars. 

The situation is somewhat better in the 7-ray band. In general, three 7-ray 
production mechanisms - the curvature radiation, synchrotron radiation, and 
inverse Compton scattering - contribute to the magnetospheric 7-radiation 
from MeV to TeV energies. Also, it is recognized that the pair cascade pro- 
cesses in the magnetosphere, play an important role in the formation of broad- 
band 7-ray spectra (Daugherty & Harding 1982). One of the key questions in 
the theory of 7-ray pulsars remains the site(s) of 7-ray production. Presently 
two type of models are discussed in this regard - the polar cap (Arons 1981, 
Daugherty & Harding 1982, Usov & Melrose 1995) and outer gap (Cheng et 
al. 1986, Romani 1996, Hirotani & Shibata 1999) models. In both models the 
primary 7-rays are produced by high energy electrons due to the curvature 
radiation. Some fraction of the curvature photons are converted into the 
pairs which initiate electromagnetic cascade in the magnetosphere. 

The winds of plasma ejected by RPP carry off the major fraction of their 
rotational energy. While initially the wind is dominated by the Poynting flux, 
at the termination shock almost all energy is believed to be in the form of 
kinetic energy of wind's bulk motion (e.g. Kennel and Coroniti 1984). The 
mechanism which could provide such an efficient transformation of the rota- 
tional energy into the kinetic energy of the cold pair wind is unknown. The 
identification of this mechanism remains a challenge for the pulsar physics 
(Michel 1982, Coroniti 1990, Lyubarskn & Kirk 2000). 

The bulk Lorentz factor of the unshocked wind is believed to be within 10^ — 
10''. Although the wind is magnetized, it does not emit synchrotron radiation 
because the electrons of the wind move together with the frozen into the 
plasma magnetic field. At the same time such a wind could be directly observed 
through its Inverse Compton (IC) radiation. In fact, the inverse Compton 
emission of the wind is unavoidable because of bulk motion Comptonization 
by external low-energy photons of different origin. It has been shown that the 
illumination of the wind's base by soft thermal photons from the surface of 
the neutron star and by nonthermal radiation from the magnetosphere should 
be accompanied by noticeable IC 7-ray emission (Bogovalov & Aharonian 
2000). The IC photons are expected in the energy interval between 10 GeV 
and 10 TeV, depending on the wind's Lorentz factor. 

The pulsar winds terminated in the interstellar medium result in strong shocks 
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which lead to the formation of synchrotron and IC nebulae around the pulsars. 
The broad-band nonthermal emission of the Crab Nebula (De Jager & Harding 
1992, Atoyan & Aharonian 1996), calculated within the model of Kennel & 
Coroniti (1984), is in general agreement with observations. Not that while 
the spectrum of IC radiation of the unshocked wind is determined primarily 
by the wind's Lorentz factor, the latter has less (direct) effect on the broad- 
band spectrum of the plerion. There are good reasons to believe that faint 
X-ray synchrotron and IC nebulae should surround many other pulsars as 
well (Aharonian et al, 1997). 

Below we discuss the expected characteristics of 7-radiation from these three 
regions for the case of Vela pulsar, the brightest GcV 7-ray source on the 
sky which can serve as a standard candle for the low threshold ground-based 
instruments in the Southern Hemisphere. 



3.1 Magnetospheric multi-GeV radiation from Vela 

All pulsar models predict that the pulsed radiation above 5 GeV is domi- 
nated by the curvature radiation of ultrarelativistic electrons accelerated in 
the quasi-static electric fields. It is generally believed that the charge den- 
sity in the pulsar magnetosphere is high enough to screen the electric field 
£"11 parallel to the magnetic field, thus the corotation condition E • B = is 
maintained everywhere except for a few locations. The regions where £"11 7^ 
can exist quite close to the surface of the neutron star - in the polar caps 
(polar cap/inner gaps model), or at a distance comparable to the light cylin- 
der along the null charge surface defined by the condition ■ B = 0, where 
the corotation charge density changes sign [outer gap model) (see, however, 
Hirotani & Shibata 2001). 

In the polar cap model (Arons 1981, Daugherty & Harding 1982) acceleration 
of primary electrons occurs at distances of an order of the neutron star ra- 
dius. These electrons are accelerated up to ~ 10^'^ eV and produce curvature 
radiation (CR) with energy of about 10 GeV. Although this radiation ini- 
tially is emitted along the magnetic field lines, during the motion of 7-rays in 
nonuniform magnetic field the angle between the wave vector of the curvature 
photons and the magnetic filed line becomes quite large, so the 7-rays inter- 
acted with the B-field to produce electron-positron pairs. The optical depth of 
the pulsar magnetosphere drastically increases with 7-ray energy. In the polar 
cap model this leads to a very sharp cutoff around 10 GeV. 

In the outer gap model (Cheng et al. 1986, Romani 1996) the electrons are 
accelerated at larger distances from the neutron star with essentially reduced 
magnetic field, therefore the pair production process becomes less significant. 
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Fig. 3. Gamma-ray spectra from Vela. The pulsed magnetospheric MeV/GeV emis- 
sion is shown by thick solid lines. The outer gap model (marked as "1") predicts 
exponential cutoff of the pulsed radiation above 10 GeV. However, in the TeV 
range this model predicts an additional pulsed IC component the intensity of which 
strongly depends on the spectrum of the magnetospheric IR emission. Solid curve 
in the 10 TeV region corresponds to the IR emission with power-law index 1.8 and 
a low energy cutoff at 0.01 eV, while the dotted line corresponds to the IR cutoff 
at 0.001 eV. The spectrum predicted in the polar cap model (marked as "2") has 
sharp super exponential cutoff below 10 GeV. For comparison the EGRET mea- 
surements are also shown - the points with thick error bars indicate the total 7-ray 
fluxes, and the fluxes of the unpulsed component are shown with thin error bars. 
The flux of tentative detection of TeV gamma-rays by the CANGAROO collabora- 
tion (Yoshikoshi et al. 1997) is also shown. The solid line is the extrapolation of the 
unpulsed EGRET spectra to TeV energies (Fierro et al. 1998). 

The high energy cutoff in the 7-ray spectrum is defined by the exponential 
cutoff of the curvature radiation and by the high energy tail of the electron 
spectrum. For the truncated electron spectrum, the resulting 7-ray spectrum 
has the following form 



where Ec is the cutoff energy of the curvature radiation. 

The predicted spectra of pulsed 7-ray emission by the polar cap and outer gap 
models (e.g. Thompson 2001) are shown in Fig. 3. The absolute fluxes in both 
cases are normalized to the fluxes reported by EGRET below several GeV 
(Kanbach et al. 1994). The corresponding differential count rates are shown 
in Fig. 4. The 7-ray detection rates are calculated from 






(12) 



14 



> 

(L> 
O 
\ 

m 

\ 

CO 

C 
=1 
O 

o 

(» 



C3 




0.01 



Fig. 4. Detection rates of magnetosheric 7-rays from Vela by 5@5 expected within 
the outer gap (curve 1) and polar cap (curve 2) models. The background count rates 
from cosmic rays calculated for 3 geomagnetic cutoff energies at 0, 10, and 20 GeV 
are also shown (dot-dashed lines). 

where M{E — E') is the energy spread function; the acceptance factor 0.67 is 
due to the collection of photons within the angle equal to the 1 sigma angular 
resolution. The background count rates are defined as 

R(E) = J An(E) ■ J(E')Aeff{E')M{E - E')dE', (13) 



where AQ{E) = 27r(l — cos^). The calculations of the integral count rates 
above 1.8 GeV yield the following values: Rpc = 21 s~^ and Rqq = 38 s~^ for 
the polar cap and outer gap models, respectively, the background detection 
rates are R^ = 5.6 and 1.3 s~^ for the 10 GeV and 20 GeV geomagnetic cutoffs, 
respectively. 

The minimum time for detection of a signal at 5a statistical significance level 
is defined from {Non — Nojf)/ ^ Non + Noj / = Vi{Rj/ ^2Rc + R.y) = 5, where 
Non is the number of counts in the direction to the source, and Noff is the 
counts in the "off direction" . According to this condition, a Scr DC-signal from 
the Vela would be detected by 5@5 after 0.9 sec and 1.8 sec observations of 
fiuxes predicted by outer gap and polar cap models, respectively, assuming 
that the observations are performed at zenith angles corresponding to approx- 
imately 10 GeV geomagnetic cutoff. The difference in the total count rates 
38 — 21 = 17 s~^ in these two models can be measured on the level of 5a after 
approximately 3 sec of observation. 

The 5@5 type instruments should detect GeV 7-rays from many other pulsars 
as well. To estimate qualitatively the minimum detectable fiuxes we assume 
that the spectrum of Vela is more or less typical for other pulsars, and accept 
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Fig. 5. Density of soft radiation components at the light cyHnder of Vela pulsar. 
Observational points of optical emission are from Nasuti et al. (1997). The thermal 
emission is approximated by blackbody spectrum with T = 9.2 • 10^ K (solid line). 
The nonthermal emission component is approximated by a power-law with photon 
index 1.8 (dashed line). 

a conservative estimate for the detection rate calculated for Vela within the 
polar cap model, Rpc = 21. Then, we arrive at a conclusion that 3 h obser- 
vations time should be sufficient to detect unpulsed (DC) 7-ray signals from 
pulsars at the level of 3 to 10 mVela depending on the available observation 
angles and corresponding geomagnetic cutoffs. The pulsars with energy spec- 
tra steeper than the spectrum of Vela would require more time to achieve such 
a sensitivity. At the same time, for pulsars with spectra harder than Vela in 
the energy region above 5 GeV, a deeper probe can be performed for the same 
observation time. 



3.2 Magnetospheric TeV gamma rays ? 

The high energy 7-rays produced in the inner magnetosphere are heavily 
absorbed before they escape the pulsar, therefore we do not generally expect 
magnetospheric TeV 7-ray emission except for millisecond pulsars that have 
relatively modest surface magnetic field of about 10^ G (see e.g. Bulik, Rudak 
& Dyks 2000). Within the outer gap models, non- negligible TeV 7-ray fluxes 
from pulsars like Vela have been predicted by Romani (1996) and Hirotani 
(2001). The basic idea is the following. The electrons responsible for the GeV 
curvature radiation should have Lorentz factor of about 10^ or more, thus in 
the case of sufficient photon fields these electrons could radiate also through 
the inverse Compton scattering. At ffist glance, although the synchrotron 
radiation of secondary electrons may serve as a natural target to support this 
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channel of 7-radiation, but it has been pointed out by Kwok et al (1991) that 
the TcV electrons have a small chance colliding with the synchrotron photons 
because they both move in the same direction. Thus, the production of TeV 
emission of inverse Compton origin requires other target photon field(s) which 
would cross the trajectory of primary electrons. Also, the magnetic field in 
the production region should be rather weak, B < 10^ G (Usov 1983) in order 
to avoid the 7-ray absorption. If these two conditions are somehow fulfilled 
in the GeV 7-ray production region, it is then easy to estimate the expected 
TeV 7-ray flux. 

The outer gap model predicts that the bulk of 7-rays is produced through 
the curvature radiation. The spectrum of this radiation is characterized by a 
high energy cutoff at = Shc-f^ /2Rc, defined by the electron Lorentz factor 
7 and the curvature radius Rc of the magnetic field line (Ochelkov & Usov 
1980). Therefore the Lorentz factor of electrons can be estimated from cutoff 
energy Ec in the observed 7-ray spectrum, provided that the cutoff does not 
have other origin, i.e., it is the result of 7-ray absorption. The intensity of 
the curvature radiation, Equ = 2ce^7^A^"c/3/?^, depends on the Lorentz fac- 
tor, the total number of accelerated electrons N^, and on the curvature radius. 
Some of the curvature photons can be converted into secondary e"^ pairs which 
in their turn would emit synchrotron photons and thus modify the initial CR 
spectrum. But the full energy flux of 7-rays, despite of possible spectral modifl- 
cations, remains close to Ecr because the total energy of emission is conserved 
(we neglect here the energy left with the cooled secondary electrons). This im- 
plies AQd^ J F{E)dE ~ -Ecr • where d is the distance to the source and AQ is 
the radiation beaming solid angle. This condition allows us to find the param- 
eter A^e from the observed 7-ray luminosity of the pulsar with an uncertainty 
contained in the product Afld^. Although the combined uncertainty in this 
parameter could be as large as a factor of 10, it does not have any effect on 
the estimate of the magnetosphcric TcV fiux if we normalize the calculations 
to the observed MeV/GeV fiux. Indeed, the same electrons responsible for the 
curvature radiation in the MeV/GeV range, emit also 7-rays through the IC 
channel boosting the soft ambient photons to TeV energies. The spectrum of 
this radiation is calculated with dN/dE^ = J a{E^,'~f){l — vc)n(e)phcA"c(ie, 
where a{E^,^) is the differential cross section of the inverse Compton scat- 
tering of an electron with Lorentz factor 7 and velocity v interacting with a 
photon with energy e and velocity of motion c; nph{e)de is the density of seed 
photons. 

While the parameter is derived from the observed MeV/GeV 7-ray en- 
ergy ffux, the electron Lorentz factor 7 is estimated from the observed cutoff 
in the 7-ray spectrum. Thus, we can write a simple relation between the re- 
sulting time-averaged energy fiux of magnetosphcric TeV 7-rays of IC origin, 
dFjc/dE^, and the observed energy flux Fqr of 7-rays integrated from 100 
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MeV to 30 GeV, 



dFic 



where Ac is the electron Compton wavelength. Apparently, the light curve 
of the TeV IC radiation follows exactly the light curve of the MeV/GeV CR 
radiation. Therefore, as mentioned above, this estimate of the IC TcV flux 
does not depend on the actual angular distribution of 7-ray emission and 
on the distance to the source if we normalize the calculations to the observed 
energy flux in MeV/ GeV 7-rays. 

At the same time, the IC photon flux significantly depends on the density 
and the spectrum of soft seed photons. The densities of these photons at the 
light cylinder of Vela are shown in Fig. 5. It consists of two components. The 
thermal emission of the neutron star is described by a blackbody spectrum. 
Several fits were proposed for this component (Ogelman et al. 1993, Page 
et al. 1996, Seward et al. 2000). Here we use the following parameters T — 
9.2 • 10^ K, = 10 km, assuming d = 250 pc. In Fig. 5 we take into account 
the possibility for the thermal photon density to fall off with distance to the 
star faster than since the angular distribution of the photons tends to be 
mono-directional at the motion from the star. This effect decreases the photon 
density at the light cylinder by a factor of 2. At larger distances the photon 
density decreases as r^^. The power-law (most likely, nonthermal) component 
of soft emission dominates over the thermal component at optical and longer 
wavelengths. Optical emission from Vela was measured by Manchester et al. 
(1978), Ogelman et al. (1989) and Nasuti et al. (1997). The densities of optical 
photons at the fight cyfinder of Vela based on the measurements of Nasuti 
et al. (1997) are shown in Fig. 5. The recent observations of Koptsevich et 
al. (2001) with VLT show that the optical emission seems to be smoothly 
continued into the near IR region down to 1 eV. Unfortunately, no data below 
1 eV is available. The e~^-^ power-law approximation to the observed optical 
data is shown in Fig. 5 by the dashed fine. 

Fluxes of the magnetospheric TeV radiation, calculated for two different as- 
sumptions for the position of the low-energy IR cutoff, are shown in Fig. 3. 
While the cutoff at 0.01 eV gives 7-ray spectrum with a distinct line- type 
feature around 20 TeV, the infrared cutoff at 0.001 eV makes the spectrum 
broader with a maximum shifted to several TeV. Although the expected fiuxes 
seem to be detectable by 5@5, the search for such a pulsed multi-TcV com- 
ponent of magnetospheric radiation could be done earlier (and perhaps more 
effectively) by the lACT arrays fike CANGAROO-3 and H.E.S.S. that are 
better designed for detection of TeV 7-rays. 



18 



3.3 Radiation from the unshocked wind 



The ultrarelativistic winds from the Crab and Vela pulsars are believed to be 
essentially cold, therefore the regions of the unshocked wind outflows should be 
under-luminous (no synchrotron radiation). This is well seen in the Chandra 
X-ray images of the nebulae surrounding the Crab (Weiskopff et al. 2000) and 
Vela (Helfand et al. 2001) pulsars. However these regions do not have to be 
similarly dark at higher frequencies. In fact, the unshocked ultrarelativistic 
winds can produce non-negligible high energy 7-radiation via IC radiation 
of ultrarelativistic (in the frame of observer) electrons moving through the 
ambient photon fields. Even in the case of isolated pulsars, an adequate photon 
density of target photons of both thermal and nonthermal origins (see Fig. 5) 
can be supplied by the pulsar itself. 

The radiation features of the unshocked wind in the Crab has been discussed 
by Bogovalov & Aharonian (2000). A similar picture should occur in the wind 
from the Vela pulsar, although the structures of the winds in these pulsars 
need not be identical. The recent observations of the compact and rather 
weak nebula surrounding the Vela pulsar indicate that the ratio of the Poynt- 
ing flux over the kinetic energy flux (the so called a magnetization parameter) 
is close to 1 (Helfand et al. 2001), while in the Crab it should be of order 
10^'^ (Kennel and Coroniti 1984). This parameter describes the ratio of the 
Poynting flux to the kinetic energy flux of the wind. The Poynting flux in 
the wind is equal to crii^kin, where E'kin is the total kinetic energy flux of the 
wind, and correspondingly ii^kin = -E'rot/(l + cr)- The longitudinal distribution 
of the wind Lorentz factor, after the acceleration, follows directly from the 
energy conservation law. In the simplest case of a wind with isotropic density, 
the Lorentz factor of the wind can be presented as T{a) = Tq + rmaxCOS^o; 
(Bogovalov 1999), where Fq ~ 100 is the initial Lorentz factor of the plasma 
produced in the pulsar magnetosphere (Daugherty & Harding, 1982), a is the 
longitude above the rotational equator of the pulsar. The recent Chandra ob- 
servations of Vela give a = 34.6° (Helfand et al. 2001). This implies that even 
for the wind with an (quasi) isotropic density the energy flux is concentrated 
at the equatorial plane which explains (Bogovalov & Khangoulian 2002) the 
observed toroidal structures in the inner parts of the Crab (Weiskopff et al. 
2000) and Vela (Helfand et al. 2001) nebulae. 

An important feature of the unshocked wind is that the electrons have, under 
certain conditions, nonzero toroidal velocity. The energy and angular momen- 
tum conservations result in the relation (Bogovalov & Tsinganos 1999) 

, rn Fo 
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Here is the angular velocity of the pulsar. It can be shown that this equation 
is valid also for the plasma flow in an arbitrary non-axisymmetric magnetic 
field. Therefore, if the plasma is somehow accelerated to a Lorentz factor 
r ^ Fq, the azimuthal angle of the wind is equal to the ratio v^p/c — Rl/t 
(Bogovalov & Aharonian 2000), whether the wind is kinetic energy dominated 
or not. Here i?L is the light cylinder, and r is the distance to the axis of 
rotation. From this equation follows that the acceleration of the wind in the 
equatorial plane to F ^ Fq is possible at distances exceeding the radius 
of the light cylinder. Also, this equation implies that there is always nonzero 
angle 9 between the velocity of electrons and the radial direction of photons. 
This angle is defined as tan^ = Rh{r'^ — Rl)^^^'^ cos a. Therefore, the IC 
scattering of electrons on the soft thermal and nonthermal emissions of the 
neutron star/pulsar inevitably results in production of high energy 7-rays. 

The spectral features of the wind emission depend on the particle flux in the 
wind N and the Lorentz factor as 

4in = A^(ro + ^F^a.)meC^ (16) 

We may expect that Fmax of the wind from the Vela pulsar lies in the same 
range as in the Crab pulsar (10^ — 10^). Here we treat F^ax as a free pa- 
rameter, and determine N through Eq. (16). For the Vela pulsar with E^^t — 
7 ■ lO'^^erg/s, the magnetization parameter was assumed a = 1. 

In Fig. 6 we show the spectra of the IC radiation of the wind for several 
values of Fmax- The IC emission of the wind consists of two "thermal" and 
"nonthermal" components produced via the inverse Compton scattering of 
the wind electrons on the thermal and nonthermal soft photons of the neu- 
tron star/pulsar, respectively. The spectra of these two components are rather 
different. Because of the scattering of the monoencrgctic wind electrons on 
the hot thermal photons with narrow Planckian distribution, the spectrum of 
the "thermal" IC component has a specific shape with a pile-up (thin solid 
curves). The latter is especially pronounced for large Lorentz factors of the 
wind because the Compton scattering occurs in deep Klein-Nishina regime. 
The broad spectrum of the nonthermal soft radiation makes the resulting IC 
7-ray spectrum also very broad (thin dashed curves). Although at low energies 
this component dominates over the "thermal" IC component, fortunately it 
cannot mask the pile-up of the "thermal" IC component that contains direct 
and unambiguous information about the wind Lorentz factor. This is explained 
by the dominance of the thermal radiation over the soft nonthermal emission 
above 10 eV (see Fig. 5). Both 7-ray components are generally expected to be 
modulated with the pulsar period. While in the case of axially uniform wind 
the periodic character of the "thermal" IC radiation could be significantly de- 
stroyed, but the "nonthermal" IC component should be always periodic. This 
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Fig. 6. The "thermal" (thin solid curves) and "nonthermal" (thin dashed curves) IC 
gamma-ray emission components from the unshocked wind of Vela. It is assumed 
that the wind with Lorentz factors Fmax = 

6-10^ (1), r,„ax = 2-105 (2), r„,ax = 6-105 
(3) and Tmax = 2 - 10^ (4) is accelerated at = i?L cos a. The heavy solid curves 
correspond to the total - "thermal" plus "nonthermal" 7-ray emission. 3a sensitivity 
threshold of 5@5 for a 3-h observation time is calculated for two (0 and 10 GcV) 
values of the geomagnetic cutoff. The sensitivity of GLAST corresponding to 30 d 
observation time is also shown. 



constraint gives us an additional possibility to separate these two components. 
And finally, it should be noted that the "nonthermal" IC component can be 
formed only if the nonthermal soft emission has magnetospheric origin. Oth- 
erwise, i.e. if the production region of the soft nonthermal photons coincides 
with the wind acceleration region (which formally cannot be excluded), the 
"nonthermal" IC component would be strongly suppressed because both the 
electrons and photons move in the same direction. 



The 7-ray fluxes of the unshocked wind shown in Fig. 6 are calculated assuming 
that the wind is accelerated at the minimum possible distance from the pulsar 
light cylinder, i?„ = Ri^/ cos a. We also show the flux sensitivities of 5@5 
corresponding to two different assumptions on the geomagnetic cutoff (0 and 
10 GeV). It is seen that in either case the sensitivity and the energy resolution 
of 5@5 should be sufficient to detect the distinct pile-up in the "thermal" 
IC radiation spectrum and to discriminate it from the broad-band emission 
components associated with the "nonthermal" IC radiation of the wind, as well 
as from the broad-band smooth IC radiation originating in the nebula (see 
Sec. 3.4), provided that the wind is formed not far from the light cylinder. 
Unfortunately, even at this most favorable condition, the detection of the 
radiation of unshocked wind of Vela seems not achievable for GLAST. 

In Fig. 7 we show the detection rates of IC 7-rays from the unshocked wind of 
the Vela pulsar calculated for two sites of the wind acceleration (more specif- 
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Fig. 7. Detection rates of the total ("thermal" plus "nonthermal") IC 7-radiation 
of the unshocked wind. The Lorentz factors of the wind are the same as in Fig. 6. 
The solid lines correspond to the position of wind acceleration at = Rh/ cos a, 
the dashed lines are calculated for = 2i?L- 

ically, in the sites of transition of the wind from the Poynting-flux-dominated 
to kinetic-energy- dominated regime) — — Ri,/ cos a (solid curves) and 
i?w = 2i?L (dashed curves). Although the difference in these two values of i?w 
is quite small, the difference in expected fluxes and corresponding detection 
rates is significant. From results presented in Figs. 6 and 7 we may conclude 
that the detection of the IC radiation from unshocked wind of the Vela pulsar 
is possible only if the wind acceleration takes rather close to the light cylinder, 
otherwise the IC 7-ray flux would be strongly suppressed. Are there special 
reasons to believe that the wind is accelerated so close to the light cylinder 
? We obviously do not know the exact answer to this question. However, we 
may refer the reader to Chiueh, Li & Begelman (1998) who concluded, from 
general analysis of an ideal MHD flow, that the acceleration of the wind to a 
low (J < 1 state is possible only in the immediate neighborhood of the hght 
cylinder. Bogovalov (2001) has argued that the wind could be accelerated to 
cr ~ 1 at the light cylinder due to the magneto-centrifugal mechanism. But 
these are theoretical arguments. What is most needed is the direct observa- 
tional evidence. The results in Figs. 6 and 7 show that 5@5 should be able to 
give an answer to this important question. 

3.4 IC gamma-radiation from the synchrotron nebula 

The termination of pulsar winds leads to randomization of velocities of rcl- 
ativistic electrons and formation of bubbles of ultrarelativistic magnetized 
plasma called pulsar driven nebulae or plerions. The nebulae surrounding the 
Crab and Vela pulsars are the most prominent representatives of plerions in 
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our Galaxy. The wind electrons injected into the Crab Nebula radiate almost 
20% of the pulsar's rotational losses in the form of synchrotron X-rays. The 
size of the X-ray nebula is constrained by the short lifetime of multi-TeV elec- 
trons caused by their synchrotron coohng. The nebula in X-rays above 0.1 keV 
has a characteristic toroidal shape with two jet-hke features along the axis of 
rotation. The linear size of the X-ray toroidal structure is about 1.2 ■ 10^^ cm. 
The basic properties of the Crab nebula are well described by the hydrody- 
namical interaction of the cold supersonic pulsar wind with the interstellar 
medium (Kennel & Coriniti 1984). 

The plerion in Vela consists of a relatively modest compact ncbiila of angu- 
lar size of 1 arcminute. For the distance to the source of about 250 pc this 
corresponds to a linear size of 2 x 10^'' cm. The recent Chandra observations 
revealed the structure of the compact nebula which appeared quite similar to 
the Crab with the same toroidal shape in the equatorial region and jet-like 
features along the pulsar's rotation axis. At the same time, the luminosity of 
the Vela nebula in the 0.1 -10 keV band is only 5 ■ 10~^ of the spin down 
luminosity of the pulsar. The small size and the low X-ray luminosity of the 
compact nebula was interpreted as indication of a large magnetization param- 
eter of the wind in Vela, cr 1 (Helfand et al. 2001), although alternative 
interpretations cannot be excluded. 

The fluxes of unpulsed 7-ray emission reported by EGRET in the energy 
interval between 100 MeV and 10 GeV from Vela are shown in Fig. 3. The 
limited angular resolution of EGRET (of about 1.5° at 1 GeV) does not allow 
any conclusion concerning the size of the 7-ray production region. Because 
of the low X-ray luminosity of the compact nebula, Fierro et al. (1998) have 
suggested that the unpulsed GeV emission is produced by the pulsar itself. 
However, this radiation could be produced in a larger region extended up 
to 1 degree or so. The Vela pulsar and the compact nebula are embedded 
in a hard X-ray (2.5 -10 keV) structure — the so-cahed Vela X. If Vela- 
X is energized by the pulsar wind, the pressure in this bubble with volume 
V = AirR^/S should be Ej-ott/V, t = 10^ yr being the age of the pulsar. For 
the pressure in the surrounding remnant of the supernova explosion Pg^t ~ 
8.5 • lO-^'^erg/cm"^, (Markwart and Ogelman 1997), the radius of the bubble 
is estimated to he R — 8.4 • 10^^ cm which corresponds to 0.64° angular 
radius, i.e. quite close to the angular size of Vela-X (Willmore et al. 1992). 
Can this interesting coincidence be interpreted as a hint that Vela-X is in fact 
a plerion harboring the energetic particles ejected by the pulsar? Interestingly, 
the power-law extrapolation of the spectrum of unpulsed GeV emission (Fierro 
et al. 1998) towards higher energies matches the tentative TeV flux reported 
by the CANGAROO collaboration (see Fig. 3). Given the large uncertainties 
both in the flux and position/size of the TeV signal, we perhaps should not 
overestimate the signiflcance such a coincidence. Moreover, the TeV flux shown 
in Fig. 3 is too high to be explained by the most natural 7-radiation mechanism 
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Fig. 8. Detection rates of the unpulsed radiation from Vela (solid line). The back- 
ground count rates of cosmic rays within Q = 0.5° calculated for 3 different values 
of the geomagnetic cutoff (0, 10, and 20 GeV) are also shown. 



in the nebula, by the inverse Compton scattering on 2.7 K CMBR, unless 
we assume unusually low magnetic field in the production region. Also, the 
energy spectrum of 7-rays predicted by this mechanism is quite different from 
a simple power-law approximation in the broad energy region. 

Nevertheless, in order to demonstrate the detection capability of 5@5 it is in- 
teresting to adopt, as a working hypothesis, that the radiation of the extended 
region in Vela is described by a power-law spectrum of 7-rays extending from 
GeV to TeV energies as shown in Fig. 3. The expected differential count rates 
corresponding to this extrapolation are shown in Fig. 8. The integral count 
rate of the unpulsed emission is estimated as 1.7 s~^. In Fig. 8 we also show 
the background rates caused by cosmic rays within the 0.5 "angular cone (the 
angular radius of Vela X) for three values of the geomagnetic cutoff in the elec- 
tron spectrum. In particular, for 10 GeV geomagnetic cutoff the corresponding 
integral background count rate is 20 s~^ which implies that the unpulsed emis- 
sion from Vela-X can be detected at the 5(7 statistical significance level for the 
observation time of about 6 min. 

The large 7-ray photon statistics and good angular resolution should allow 
spectroscopic measurements and determination of the size of extended pleri- 
onic component of 7-radiation in a broad energy region extending from several 
GeV to ~ 1 TeV. Note, how^cvcr, that such important spectrometric and mor- 
phological studies above 100 GeV can be performed already in the near future 
by the H.E.S.S. and CANGAROO-3 lACT arrays. At lower energies down to 
1 GeV, GLAST will contribute key information. 
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4 Search of periodical signal from radio-quiet RPP 

The nature of 170 unidentified point-like sources of > 100 MeV 7-rays detected 
by EGRET (Hartman ct al. 1999) remains mostly unknown. It is likely that a 
significant fraction of these objects belongs to RPPs with radio beams which 
do not cross the Earth orbit (Grenier 2000). The beams of 7-ray emission 
are believed to be much wider than the radio beams (see e.g. Romani and 
Yadogaroglu, 1995; Harding and Zhang, 2001). Therefore, a chance to detect 
7-rays from a pulsar can exceed the chance to detect radio emission. Obviously 
the most direct way to identify the unidentified EGRET sources with pulsars 
would be the detection of pulsed 7-radiation. 

The search for periodic signals in the 7-ray data is not a simple issue (see e.g. 
Chandler et al. 2001). The pulsed components of a// EGRET pulsars have been 
discovered only after applications of information about the period and phase 
obtained at other energy bands. However, in many cases the time structure 
at different wavelengths could be different. Therefore, it is crucial to search 
for periodic 7-ray signals without relying on observations at other energy 
bands. In turn, this search requires an adequate 7-ray photon statistics. The 
statistics of photons increases with observational time, and in principle could 
be very large, in particular for large-field-of-view instruments like GLAST. 
However, the valuable observation time can be increased infinitely only if the 
source of the periodic signal is a perfect clock. Pulsars are not perfect clocks. 
In fact, their periods increase with time due to the rotational losses. The 
phase of the signal varies with time as (/? = 27r(tz/ + t^u), where u = 1/P 
and P is the period of rotation, t is the time of observations. The variation of 
the phase due to the change of the period of rotation can be neglected while 
t^v < 1/2. Thus, the duration of the total time interval for the search of a 
periodical signal in absence of an a priory information about the period and 
period derivative should satisfy the condition Tobs < ^l/2z> . For example, the 
first derivative of the rotational frequency of the Crab pulsar z> = 3 ■ 10^^'^ s^^ 
requires observation time Tq^s < 4 • 10*^ s ~ 10 h. On the other hand, the 
number of detected photons should be large enough to allow determination of 
the period. Let's assume that the fight curve of the source consists of a single 
narrow pulse. This assumption agrees with the observed pulse shapes above 5 
GeV for Crab, Vela, PSR 1706-44 and Geminga pulsars (Thompson 2001). In 
this case the epoch folding method allows accumulation of i?(> 5 GeV)Tobs 
pulsed photons. Let's assume that the full phase consists of uq bins, and that 
the pulse is located within one of these bins. Then, the m-sigma detection of 
the 7-ray signal in one trial requires i?min > ?7T.y'i?b/^oT'obs, where R\, is the 
total background count rate, Uq is the number of bins in the light curve. Below 
we adopt 5(7 condition for detection of the periodical signal. This implies that 
the probability to detect signal by chance multiplied to the number of the 
trial fight curves should be less then 3 • 10~^. The number of the trial light 
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curves A^tr = (j^max ~ ^min)/Si', wheie i^max = 10^ Hz is the maximum possible 
frequency (corresponding to a millisecond pulsar) and z/„it„ = 1 Hz is the 
minimum possible frequency, = l/Tghg is the frequency separation. Thus, 
A^tr ~ 10^ for the accepted parameters and therefore the probability to detect 
the pulse by chance in one trial should be less then 3 • 10~^^. This means that 
the 7-ray signal should be detected at the level of m—7.5. 

It is believed that GLAST will be able to detect as many as 250 pulsars, half of 
which would be unknown at other wavelengths (McLaughlin & Cordes 2000). 
Also, GLAST can reveal the pulsed emission component from a number of 
unidentified EGRET sources if they have indeed pulsar origin. The same can 
be done by 5@5 type instruments. Moreover, 5@5 should be able to identify 
the future "unidentified GLAST" sources with pulsars if the spectra of these 
objects extend beyond 5 GeV. For the background count rate Ry, — Re — 
5.6 s~^, and assuming Tq^s = 10^ s ~ 3 h (a standard time of observations of 
a source per night) and Uq = 10, for m = 7.5 we obtain the 7-ray count rate 
-Rmin = 0.06s^^. This corresponds to the fiux level of approximately 3 mVela if 
we adopt the minimum possible count rate R ~ 20s~^ from Vela as predicted 
by the polar cap model. 



5 Irreguleir Veiriability of the 7-ray radiation from radio pulsEirs 

The 5@5 telescope could be quite effective for studies of 7-ray flux variability 
from radio pulsars. Presently at least two types of irregular time behavior of 
radio pulsars are established. One of them is the so-called "timing noise" in 
which the pulse phase and/or frequency of radio pulses drift stochastically 
(Helfand et al. 1980, Cordes & Helfand 1980). Another type of irregularity in 
the radio emission from radio pulsars are "ghtches" - sudden increase in the 
pulse frequency produced by apparent change in the momentum of inertia of 
neutron stars. These rather rare events occur mostly in young pulsars (Reichly 

6 Downs 1971, Lyne 1996). Ghtches were observed from 27 pulsars (Wang et 
al. 2000). 

The Vela pulsar seems to be the best object to study the irregular variability 
of pulsars in 7-rays by 5@5. A dozen glitches have been recorded over the 
past three decades from this source. EGRET data show that the fiux in the 
range above 100 MeV is variable. Remarkably, there is evidence of variation 
of the spectral index of the radiation, as well as of the relative intensity of two 
peaks of the light curve (Ramanamurthy et al. 1995). The variability of 7-ray 
emission from the Vela pulsar detected by EGRET has been derived for 14-day 
observations but the photon fiux was too low to investigate the variability on 
smaller time scales. 5@5 is expected to allow detection and resolution of the 
7-ray variability on sub-hourly timescales. 
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6 Summciry 



The high flux sensitivity of 5@5, supported with unprecedented photon statis- 
tics in one of the most informative windows of electromagnetic radiation for 
the pulsar physics, E > 5 GeV, should allow detailed studies of spectral and 
temporal characteristics of the 7-ray pulsars detected by EGRET, as well 
as very effective searches for pulsed components from unidentified EGRET 
sources. The energy and angular resolutions of 5@5 can provide unique stud- 
ies of other radiation components of pulsars associated with the unshocked 
winds and the synchrotron nebulae (plerions) in a broad energy band extend- 
ing from several GeV to 1 TeV. We believe that the new generation sub-10 
GeV ground-based 7-ray instruments like 5@5 can provide crucial insight 
into the physics of pulsars and their interactions with surrounding interstellar 
medium. 
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